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Abstract The non-perturbative domain of QCD contains confinement, chiral sym- 
metry breaking, and the bound state spectrum. For the calculation of the latter, the 
Coulomb gauge is particularly well-suited. Access to these non-perturbative properties 
should be possible by means of the Green's functions. However, Coulomb gauge is also 
very involved, and thus hard to tackle. We introduce a novel BRST-type operator r, 
and show that the left-hand side of Gauss' law is r-exact. 

We investigate a possible truncation scheme of the Dyson-Schwinger equations in 
first-order formalism for the propagators based on an instantaneous approximation. We 
demonstrate that this is insufficient to obtain solutions with the expected property of 
a linear-rising Coulomb potential. We also show systematically that a class of possible 
vertex dressings does not change this result. 
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1 Introduction 

QCD in Coulomb gauge has a long history as a gauge being particularly useful for 
the calculation of the bound-state spectrum, one of the most intricate properties of 
QCD. Furthermore, a well-developed confinement scenario exists in the Coulomb gauge 
that was developed by Gribov and elaborated by one of us [U[2]- Recently i h has 
also been shown that Coulomb gauge furnishes an upper limit for the physical string- 
tension and static quark-anti-quark potential already at the level of the two-point 
functions 3J 415]. Coulomb gauge can also be understood as a gauge which in a certain 
sense is physical: Unlike e.g. in covariant gauges, perturbatively unphysical gauge-fixing 
degrees of freedom do not propagate. Finally, its structure makes it especially suited for 
calculations at finite temperature. The results obtained in this gauge can also provide 
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input to model calculations (see e. g. [6]) in the form of potentials, and is by this 
relevant to the work (e. g. [7]) of the jubilee to whom this volume and this work is 
dedicated. 

These are particular advantages of the Coulomb gauge compared to, e.g., covariant 
gauges. However, these properties do come at the price of loosing explicit covariance, 
and therefore introducing a much more complicated structure in the mathematical 
description. As a consequence, the question of renormalizability [2ll8ll9lfT0] of this gauge 
is a very hard problem. Furthermore, the lack of covariance makes explicit calculations 
cumbersome even in perturbation theory [9lll0j. 




Therefore, the determination of the Green's functions in the perturbative and non- 
perturbative domain is far less developed than in covariant gauges. Nonetheless, we 
will present an attempt at such a calculation in the far infrared regions by means of 
a method which has proven to be rather useful in Landau gauge [Hlll2lfT3lll4j . an 
asymptotic analysis in terms of critical Green's functions. We attempt no full solution 
of the infinite tower of Green's functions, as is possible in Landau gauge |13II14] , and 
will also restrict to the case of Yang-Mills theory, thus neglecting quarks. For this 
purpose, the first-order formulation will be used. A possible translation to the more 
widespread used second-order formalism can, in principle, be made using the translation 
prescription discussed in [1 5J . 

It is a remarkable fact that such calculations in Coulomb gauge in two dimensions 
can be done analytically [16]. The results of these calculations confirm the scenario, 
which will be used throughout this work. This result is, of course, due to the strong 
constraints the Coulomb gauge imposes specifically in two dimensions. In higher dimen- 
sions, the situation is much more complex, and many attempts to obtain the Green's 
functions [T7] , in particular using variational principles [18] and lattice calculations |19l 
[20ll21H22| . have been performed. 

We proceed as follows. General aspects of the Coulomb gauge and the first-order 
formalism are discussed in section and the corresponding DSEs in section (|3]). In 
section (2) we also introduce a novel on-shell BRST-type operator r and show that 
the left-hand side of Gauss' law is r-exact. The truncation employed here, described in 
section ((4]), provide a consistent result at the level of power-counting. Unfortunately, 
if a linear rising Coulomb potential is required, this solution ceases to exist, as is 
discussed in detail in section (|4.4[) . This persists even when including consistent vertex 
modifications, as is discussed in section ([5j|. The results provide an example where a 
solution exists by power counting, but is not a solution of the system. Of course, this 
may be an artifact of truncation. The results of the investigation are then summarized 
in section (|6]). 



2 General aspects of Coulomb gauge in the first order formalism 

2.1 Phase-space action 

The local Euclidean action in Coulomb gauge in phase-space or first-order formalism 
in s spatial dimensions is 
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where 

Bf = B?(A) = e ijk (djA k + X -gA, x A k ) 
is the color-magnetic field, and 

(A x B) a = f abc A b B c . 

Here c and c are the ghost pair, and D.; = Di(A) is the gauge covariant derivative, 
DjC = diC + gAi x c. All fields are understood to carry a color index in the adjoint 
representation that is contracted. The field 6 is a Lagrange multiplier that enforces the 
Coulomb gauge condition diAi = 0. Here m is an auxiliary field that represents an 
independent color-electric field. 

The Euclidean phase-space measure is 

j dAidAo dn dcdcdbexip(—S). 

If one integrates out 7r.j, one gets the second-order or Lagrangian Euclidean action in 
Coulomb gauge. 

As in the Landau gauge [23], it will be assumed henceforth that a non-perturbative 
realization of the Coulomb gauge exists which furnishes a globally well-defined and 
unbroken BRST symmetrjQ. If this is the case, an analogue to the classical Gauss' law 
can be formulated using the BRST symmetry, as explained below. As a consequence, 
certain relations between the renormalization constants of Coulomb gauge have to hold, 
in particular Zj± a — Z c , relating the ghost renormalization with that of the temporal 
gluon propagator [2J. 

2.2 On-shell action and on-shell BRST operator 

We integrate out the auxiliary &-field, so that the gauge condition diAi = is satisfied 
identically (on-shell). Writing A^ = Aj , the on-shell action 

Sos = J d s+1 x\i-K t {d Aj - D Z A ) + i(7r 2 + B 2 ) + diWic] , 

is obtained, where B = B(A T ), and D = D(A T ). 

Since the gauge condition is satisfied identically, it is appropriate to restrict the 
integration over A T to the fundamental modular region A, a region free of Gribov 
copies. The other variables are unrestricted, and the path integral is given by 

J dA T dTTdA dcdcexp(~S O s)- (1) 

In the phase-space action, Aq appears as a Lagrange multiplier, and if one integrates 
out Aq, one obtains 5(Dj7Ti), which imposes Gauss' law0 -D,;7Tj = 0. Note that the 

1 The problem of the existence of a globally well-defined BRST charge in this gauge is not yet 
solved, and similar considerations as in Landau gauge apply |23I124| . However, a construction 
as proposed in Landau gauge |24II25| may also be possible in Coulomb gauge. Furthermore, 
the existence of a conserved and vanishing global color charge |26j may also be helpful in this 
respect. 

2 In the presence of quarks one would have instead 8(DiTTi — p qu ). 
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integrand vanishes on the boundary of the Gribov region, after integrating out the 
7r-field, because the Faddeev-Popov determinant vanishes there, and hence the Dyson- 
Schwingcr equations (DSEs) are unchanged by a cut-off there |27| . 

For purposes of deriving the DSEs, it is preferable not to integrate out Aq because 
that would produce a non-polynomial action. Instead we exhibit an on-shell BRST-type 
operator which encodes Gauss' law, and provides a definition of physical observableo 
The on-shell BRST-type operator r is defined by 

rc = iDj7Tj; rAj = rin — 0; (2) 

rA = d i D l c; rc =0. (3) 

An operator similar to r was introduced in eq. (143) of [8]. These equations imply that 
r is nilpotent, r 2 = 0. Moreover it is a symmetry of Sos, 

rSos = 0, 

as one sees by writing the on-shell action in the form 

S os = J d 4 x [iir t d Aj + i(vr 2 + B 2 ) + r(cAo)] . 

The r-symmetry is compatible with the restriction to the fundamental modular region A 
of the integral |l]) over A T because rA T = 0. 

We identify physical observables with the cohomology of r. With this definition of 
observables, one sees from ((2j and (0 that the variables Aj and 7Tj are physical, while 
(Aq, c, c, Di-Ki) form an unphysical quartet. Any r-exact quantity rX has vanishing 
expectation value (rX) = 0, and is physically equivalent to 0, rX ~ 0. Gauss' law is 
encoded in the r-symmetry because, according to ([2]), Dj7Tj is an r-exact quantity, and 
is thus equivalent to 0, D^i ~ 0. The r-symmetry is not as powerful as the BRST 
symmetry because it does not determine the part of the action that is not r-exact. 
However it does imply the cancellation of the instantaneous ghost and scalar loops 
which contain the notorious energy divergences of the Coulomb gauge. 

In the BRST formalism Aj and 7Tj are not physical. However here they can be 
because the gauge is fixed on-shell. To make contact with the usual Coulomb gauge, 
in which only the transverse degrees of freedom are considered physical, we separate 
the 7r-field into transverse and longitudinal parts, 7Tj = nj — di<j), where 4> is the color- 
Coulomb potential, and diiif = 0. These fields satisfy 

rAj = nrf = (4) 
r(M _1 c) =i(<t>- M~ l p); r<j> = (5) 

where 

p = g tyJ x Aj (6) 
is the color-charge density of the transverse degrees of freedom, and we have used 

r(M^c) = M-hD m = iM-^-Didri + gAf x irf). (7) 

3 We note that this operator is affected by the same lack of knowledge concerning its global 
definition as in case of the BRST symmetry. Once more, we assume that it is globally well- 
defined. 
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Here 

M = -d i D i (A r ) = -D i {A T )d i (8) 

is the Faddeev-Popov operator, and the last equality holds because Aj is transverse. We 
see by the preceding definition that Aj and tvJ are physical, as in the usual Coulomb 
gauge, whereas ci is physically equivalent to M~ 1 p, 

cj> ~ M~ x p, 

because they differ by an r-exact term, and moreover M p is a function of the trans- 
verse degrees of freedom Aj and nj. In the usual Coulomb gauge, the color-Coulomb 
field has the value cf> — M~ 1 p. 

3 Dyson-Schwinger equations 

With the separation m = tyJ — di<f>, the on-shell action becomes 

Sos = J d s+1 x^inJd Aj + 7}[(k T ) 2 

+(cW 2 + B 2 ] + id^D.Ao - igirf(Ai x A Q ) + d l cD l <^. (9) 
Transverse propagators are expressed in terms of scalar quantities by 
(AJ"( X )AJ\0)} = 5«»J d^ k ^Ap ijik)D l A{k) 

etc., where Pij(k) = S.^ ~ kikj/\s 2 is the projector onto the transverse subspace. The 
derivation of the DSEs from the action Sos is straightforward, keeping in mind that 
the propagators are 2 by 2 matrices, 

D T AA {k) D T A „{k)\ (D AoAo (k) D Ao<p (k)\ 

Dj A (k) D^(k) J ' V D Mo (k) D^(k) J ' 

that are inverse to 

rJ A (k) rjjk)\ fr AoAo (k) r Ao0 (fc)\ 
rJ A (k) r^{k) ) ' V r Mo (k) r^(k) ) ' 

respectively. To define the line-styles for a pictorial representation, the tree-level prop- 
agators are shown in figure [T] The tree-level vertices are given in table [T] 



Table 1 The tree-level 3-point vertices. 4-point vertices will be neglected throughout. Mo- 
mentum conservation is left implicit. Fields denote the involved fields for the vertex to fix the 
association of indices, momenta, and field-type. 



Fields 


Notation 


Tree-level value 






-igf abc ((g - r) t S jk + (r- p) 3 5 ik + (p - q)^)) 


r{p)A^ b (q)Af ) (r) 


r? ALA ° abc (p,q,r) 


9f abc Pi 


TTi a (p)A* b {q)A° {r) 


7-t7r A 1 An. abc i \ 

r tj (p, q, r) 


~igf abc S,.j 


c a {q)c h {r)A^{p) 


rf A ' abc (p,q,r) 


-igf abc n 
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<^a A Tb> 8 ab p (8 ir P iP] /p 2 )/p2 

<4>'<> V\Ar -i§ a V 

<c a c"> 5 ab /p 2 

Fig. 1 The tree-level propagators of all elementary and mixed fields. 

A complete pictorial representation for the Dyson-Schwinger equations, truncated 
at one-loop level, is given in figures [2] (^-propagator, 7r-propagator, transverse gluon 
propagator, and the time- like gluon propagator) and[3] (mixed 7r-transverse gluon prop- 
agator, mixed <f>-Ao propagator, and ghost propagator). The full expressions will be 
given only for those contributions present after truncation. A full set of DSEs could 
be generated, e. g., using the method described in [29] . An independent derivation of 
these equations can also be found in [10] . No consequences of additional mixed propa- 
gators have been included, because even if they are non-zero, they do not survive the 
proposed truncation. Also, the following expressions for the DSEs will only contain the 
parts which will be needed explicitly. The justification for the truncation scheme will 
be given below. 

The DSE for rJ A reads 

r AAij( k ) = S<lb ( k2< % - kkjj 

+i9fadc 1 0^0^ ( fc - P n; ATA ° feb (P, k - Pl -k) 

+9fadCh I ^ D t4p) D t4 k ~ p)rf AoATfeb ( P , k - P , -k) 
+igf adc k J j^d c J (p)D d c l(k - P )rf cAT f eb (p, -k, k- P ). (io) 

and with tree-level vertices 

Pij(k) rJ A (k) = Py(fc) k 2 + J?Sp a (k) J d s+1 p[p lm (p) D^(p) D AaAo {k-p) 
+ Pl D Aa4> (p) D Aotj) (p - k) (p- k) m 

+Pl Dcc{p) Dcc(p — k) (p — k)m] P m j 

N denotes the adjoint Casimir of the gauge group. Upon taking the trace on the spatial 
indices one obtains the DSE for scalar quantities, 

rUk) =k 2 + ^ T / ^i p [ '-a + (M a dUp) D AoAo (k- P ) (ii) 

+ 1 ~ ( ^ 1 P)2 P 2 DAoM D Ao< p(P -k)+ 1 ~ s i *'* )2 P 2 Dcc(p) D c£ (p - k)] , 



<d> a q> b > ■ o 

<a>™> V.Q.Q.P.P.Q.P., ^fP-pyVp 

<A a A b > VWWW 5 ab /p 2 
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Fig. 2 The DSEs for the 0-propagator, 7r-propagator, the transverse gluon propagator, and 
the time-like gluon propagator. Only black diagrams are kept in the truncation. Red diagrams 
vanish for tree-level vertices. Blue diagrams do not vanish when inserting tree-level vertices, 
but are discarded in the truncation. Lines with a dot are full propagators, and open circles are 
full vertices. 



s 





Fig. 3 The DSE for the mixed 7r-transverse-gluon propagator, the mixed </>-j4o propagator, 
and the ghost propagator. Only black diagrams are kept in the truncation. Red diagrams 
vanish for tree-level vertices. Lines with a dot are full propagators, and open circles arc full 
vertices. 



where k and p denote the corresponding unit vectors. Similarly the equations for _T W7r 
and rJ A read 

r$? j (k) = 5 ab (5 ij k 2 + k i k j ) (12) 
+l9fadC J d i^ DT A C Liv)D% Ao (k - P )r^ ATA ^ b (-k, P ,k-p) 

or with tree-level vertices 

r ™ {k) = 1 + J dS+1 P S ~ 2 ^}\' P)2 Dl A {p) D AoAo (k-p) (13) 

and 

rJ A (k) = -*o + j^f+f J d s+1 P s - 2 +_^P) 2 D Av (p) D AoAo (k-p). (14) 
For r AoAo we obtain 

r AQ A {k) = j^-^ J d s+1 p [k 2 [l - [k-pf] D T AA {p) D^(p - k)] . (15) 

The equations (|15|) and (|14[) will not be needed in a form with full vertices, and thus 
these expressions are skipped. For r^j, we have 

^(i).A 2 + s A I ^Dl% j (k-p)D d / oAo ( P )rf AoATb9f (-k,p,k-p). 

(16) 
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or with tree-level vertices 

r^{k) = k 2 + | d s+1 P [i - (fc-p) 2 ]J)L(p) d AoAo ( p - k) , 

and for r$ Ao , 

rf Ao {k)=i5 ab ^+gf dae k i J ^^Dl% j (k-p)D%^k)rf A » ATb ^(-k, P ,k- P ) 
or, with tree-level vertices introduced, 

r Mo (k)=ik 2 + ^^ I J d s+1 p[l-{k-pf]Dl A {p) D Ao<p (p-k) . (17) 
Finally for the ghost we obtain 

r^^A 2 -;/^, J ^^ D l%(k-p)D d c !(p)r^ Tb9f (-k,p,k- P ), (is) 

which takes for tree-level vertices the form 

r- cc {k) = k 2 - Bj^. J d s+1 P [i - (* . p) 2 ]z?I a (p) Dco(p - fc) , (19) 

where T ec (fc) = D^(k). 

4 Truncation with tree-level vertices 

Many of the following arguments, and of the truncation, still hold also for the ansatze 
for the vertices which will be made in section [5] For simplicity, here the tree-level 
expressions will be used, and only commented on, if some argument no longer holds 
when using non-tree-level vertices. 

4.1 Instantaneous parts of propagators 

In general the propagators depend on two variables, D(k) = Z)(|k| , fco)- However in 
Coulomb gauge the propagators of the scalar degrees of freedom contain instantaneous 
parts Dj, proportional to S(t), and these, in momentum space, depend on a single 
variable D/(|k|). Here we shall identify the instantaneous parts of propagators. 

In the phase space action (J5J), the only time derivative occurs in the term iirj doAj '. 
Accordingly we suppose that the propagators of the transverse degrees of freedom con- 
tain no instantaneous parts, but that the scalar fields may have instantaneous parts 
proportional to S(t). It can be shown by appropriately integrating out fields that prop- 
agators of the scalar fields may be expressed as 

5 ab D C c{x — y) = ( (M~ 1 ) ab (x,y) ) 
5 ab iD Ao4> (x-y) = < (M~ 1 ) ab (x,y) - (Kp) a (x) (M"V) & (y) ) 
S ab D AoAo (x-y) = < K ab {x,y) - (Kp) a (x) (Kp) b (y) ) 
5 ab D^(x-y) = < (M~ 1 p) tt (x) {M- l p) b {y) >, 
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where K is the operator with kernel 

K a \x,y) = [A/- 1 (~V 2 )Af- 1 ] a "(^,y)- (20) 

We have p — givj x Aj and M — —d^D^A) so the quantities on the right depend 
only on the transverse degrees of freedom. The propagators of the scalar fields contain 
instantaneous parts because the kernel M^ 1 (x,y) = M~ 1 (x.,y)5(xQ — yo) is instanta- 
neous, and consequently so is K(x,y) = K(x,y)S(xQ — yo). 

Assuming that the transverse fields do not contribute to the instantaneous parts, 
the last equations give the decomposition of the propagators into instantaneous (I) and 
non-instantaneous (N) parts, which reads in momentum space 

Dcc(k) = D IcS (\k\) 
D Ao <t>( k ) = D IAotj) (\k\) + D NAo< p(k) 

DAaA (k) = D IAqAo (M) + D NAoAo {k) 

D <l><t>(k) = Dm^U). 

Moreover we have the equality of the instantaneous Bose and Fermi propagators 

^ o0 (]k|) = D Jcg (|k|). (21) 

which may be shown to be a consequence of the r-symmetry. 

The ghost propagator D c e(|k|) is purely instantaneous. The integral over po m (|19|l 
yields 

^- [ d P0 D^ A (p) = D =AA (\p\), (22) 



2?r „ 

where D— AA (\]i\) is the equal-time A^-A^ propagator in momentum space. The DSE 
of the ghost, 

D-^lkl^k 2 -^^ J d s p[l-(k-p) 2 ]D =AA (\p\) Uccdp-kl), (23) 

thus involves unknown functions of only one variable. 

Later on, in section [5] we will presume that the momentum dependence of the 
vertex dressing functions appearing also are purely instantaneous, i. e., only depend 
on k and not on /cq. Therefore, everything in this subsection also applies to this case. 



4.2 Instantaneous approximation and infrared limit 

We wish to explore the hypothesis that there exists an asymptotic infrared limit of the 
DSEs which is dominated by loops containing an instantaneous propagator. As Ansatz 
we suppose that we may neglect loops that contain no instantaneous propagators on 
the RHS of the DSEs for scalar propagators, eqs. (|15p - (|17|l . This is equivalent to the 
substitutions on the RHS 



Aw,(|k|,fco) 
D AoAo (\k\,k ) 
iD Ao4> (\k\,k ) 



(24) 
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Equations JT5} - CE} then yield on the LHS 
Fa a = 

— *-Tz>A = -Tec 



„2,,2 

(2 

which are all consistent with our Ansatz. We have 



k2 + ^^ / d s p[l~(k-pf}D =AA (\p\) D IAoAa (\p-k\), (25) 



which gives 



^W|k|) = -i^ (26) 



- ^ e( | k |) • w 



This implies 

D 7 ^ oAo (|k|) o iV g 2 k 2 f M , £ ., 2 



^c 2 E (|k|) 



k2 + ^r^ / ^P^-^-P^^aQpD D IAoAo (\p-k\). (28) 



The pair of equations ((23} and ([28} determine the two propagators D c g(jkj) and 
Dj AqAo (|k|), provided that D = ^4(|k|) is known. 

We now consider the DSEs for the transverse propagators. We make the same 
approximation D AqAo (k) — > Dj AoAg (|k|) into the RHS of (fT4)l which yields the integral 



/ 



dpo £>L(p) = 0, (29) 



which vanishes in the present approximation because D\^{pq, |p|) is odd in pq, as one 
sees from the tree-level term in (|14|l . and as will be verified below. This gives 

r^flki) = -fe . (30) 



We make the same substitution D AqAo — > Dj AqAq in (|13[) and obtain 

r^(|k|) = i + ^ | rfV 8 ~ 2 s + _ ( i ' ^ D =AA (\p\) D IAoAo (|k - p|). (31) 

There is no fco dependence on the RHS, so, as indicated, J^dkl) depends on the one 
variable |k|. 

There remains to consider the DSE for r AA . Upon making the substitution iD Ao( p — > 
Dec into the RHS of (|11[) in the loops that contain instantaneous pieces, one finds 
that the loop of scalar bosons Dj Ao ^(p)Di Aol p(p — k) exactly cancels the ghost loop 
Dcc(p)D C c(p — k). This is a crucial cancellation, and it holds in this truncation only 
with tree-level vertices. Beyond such vertices, this is non-trivial, and implies constraints 
which will be discussed in section [S] Of course, in the full solution this would be guar- 
anteed by the r symmetry. 

This leaves only one loop in the equation, and we obtain 

rL(|k|) =k 2 + ^ J d s p S ~ 2 s + _ ( f^ )2 -P=^(|Pl) D IAoAo (\k-p\). (32) 
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Again there is no &o dependence on the RHS, so _T^(|k|) is a function of |k| only. 

It should be noted that there is a certain ambiguity involved here: In the DSE 
(fTTj) for rJ A there is a cross term D j Ao ^{p)D N Ao ^{p — k) containing an instantaneous 
propagator (I) multiplying a non-instantaneous propagator (N), which in principle 
should be kept. However, if we keep the one instantaneous part on the RHS of the DSE 
(|17[) for r^ Ao > then r^ Ao comes out purely instantaneous, so within our approximation 
Dn4>A = 0, and we have no way to calculate a non-zero value for this quantity. This 
is self-consistent. However, this assumption may be one of the reasons for the failure 
of finding a solution with the desired properties in this truncation. 

Using the results of the last three equations, we see that the transverse _T-matrix 
is given by 

'rJ A {\k\) k 

The elements of the inverse matrix are thus 

ol A (fe , |k|) = rX(|k|)/Q 
£>L(*o,|k|) = -ko/Q 

|k|) = rJ A (M)/Q, (33) 



where 



k 2 + rJ A (\k\) r^(|k|). (34) 



We see that D Av (ko, |k|) is odd in fcrj, as asserted. 

It follows that the equal-time propagators that appear on the RHS of the DS 
equations are given by 

If t 

D=aa(\p\) = J dp D AA (p) 



and similarly 



This implies that 



J_ f . rU\p\) _ ir r^dpD ii/2 

2W po P l + rJ A (\ P \)rJA\p\) ^rJ A (\ P \)\ ' 

D=™(\p\) = ^Z^(IpI), (36) 



so there is only one unknown function of one variable. Upon taking the ratio of (|31l) 
and (|32[) we obtain 

1Di ^ kn - few - f < 37 > 



where X and Y are the RHS of (|31[) and p2[) respectively. This provides the desired 
equation for D =AA (\]i\). However we expect that Dj AoAo (|p — k|) is highly singular 
at k = p, for example like l/|p — k| 4 , so there are infrared divergences in X and Y. 
However, these cancel by rearranging the last equation to read 

4k 2 D =AA (\k\)-DZ\ A (\k\) (38) 

Ng 2 f s-2 + (p-k) 2 ,s\ D=AA(\p\) D=AA(M) 1 

In all we have three equations, ([23} and (|28[) and (|38|) for the three functions D c e(|p|), 
Dia A (\p\) and D =AA (\p\). 
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4.3 Infrared asymptotic solution 

We suppose that these three propagators have infrared asymptotic limits that are 
described by power laws, 



9 D IA(>Ao (\p\) = 



(p2)l+K V 
(P 2 ) 1 



9 2 D= AA (\p\) = 7-2^7. (39) 



4.3.1 Solution of DSE for D c 



As discussed elsewhere (2][T] we choose a solution for which the ghost propagator 
Z) C c(|k|) is more singular than the free propagator 1/k 2 at k = 0, implementing the 
horizon condition. We therefore require that the term quadratic in k on the r.h.s. of 
(|23[) is canceled, which gives 



D ^(^ = ^0T J d s p[l-Ck-p) 2 ]D =AA (\p\)[D c ,(\p\)-D c5 (\p-k\)]. (40) 

This choice of a cancellation mechanism makes the equation invariant under the renor- 
malization group (RG), as it must be. 

By comparing powers of momenta on the left and right one finds from the power 
laws (|39p that the critical exponents satisfy 

k t + 2k a = (s - 4)/2. (41) 

The integral is ultraviolet convergent provided that s/2 < 3 + + k A , and infrared 
convergent provided that s/2 > 2 + + k a , which gives 

< ka < 1- (42) 
The integral on the r.h.s. of (|40)l was evaluated in eq. (A. 17) of [27], with the resul10 

f66 a r i = N(8-i)r(i- KA )r(2 KA ) r(s/2-n A ) 

{TA> (4tt) s / 2 r(i + K A ) r(s/2-2 KA ) r(i + s/2 + KA ) ( ' 

where we have used it/ sin(TTK A ) — —r(—K A )r(l + K A ). 



4.3.2 Solution of DSE for D IAoAo 

We next substitute the power laws (|39|l into (|28p . The integral is ultraviolet convergent 
provided that s/2 < 2 + kj- + Ky, and is infrared convergent provided that s/2 > ny. 
This gives 

2k A < k v < s/2, (44) 



4 Due to the structural equivalence to Landau gauge it seems likely that a more general 
analysis, like in the Landau gauge case 1281 . is possible. 
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where we have again used (gTJ. The l.h.s. of (J28J) is given by &yb^ 2 (k 2 ) 1+2Kzl Kv . By 
the last inequality, this dominates k 2 at low k, and we conclude that the tree-level 
term in (|28p is negligible in the infrared, and (|28|) simplifies to 



D c 2 £ (|k|) " (2nY 



J d s p [1 - (k ■ pf]D =AA (\p\) D /AoAo (|p-k|). (45) 



This equation is also invariant under RG transformations. 

We substitute the power laws (|39|l into this equation, and upon equating powers 
of momenta, we again obtain (|41|l . It is straightforward to evaluate this integral which 
yields, 

(b h 2 r 1 _ N(s - 1) £(a/2 - ny) T(l + 2k a ) r(K V - 2k a ) 

{TA> 2(4tt) s / 2 F{s/2-2k A ) r(l + K V ) r(l + s/2- k v + 2k a )' 1 1 

where we used (|41|) to eliminate k^. Upon comparing with (|43|l . we obtain a relation 
between k a and Ky, 

K A r(s/2-K V )r(K V -2 KA ) _ r(i-K A ) r(s/2-n A ) 



r(i + k v ) r(i + s/2 - k v + 2k A ) r(i + k a ) r(i + s/2 + 

^.3.3 Solution of the DSE for D*L AA 

Finally we consider eq. (|38p for the transverse propagator. We shall show that at small 
k the loop integral dominates the tree-level terms. The tree-level terms on the l.h.s. 
are of order |kp 2KT = |k| 4K ^ +4 - s and |k| 2+2KT = |k| s - 2 - 4K ^ respectively. The ratio 
of these terms is |j 4 | 8ft -^+ 6_2s i From the inequality (|42p we have k a > 0, so this power 
is positive for s < 3. Thus, at small momentum the second tree-level term |k| s-2 ~ 4Kzl 
dominates the first. Moreover the dimension of the r.h.s. is |k| s_2_2Kv ', so it, in turn 
dominates the second tree-level term by the inequality (|44[l. Ky > 2k A . We conclude 
that the loop integral in (|38|) dominates both tree-level terms in the infrared asymptotic 
limit, and this equation simplifies to 

o - $ / *v I. - > + V - m D IM . ,| k - Pi) [Ml - Ml] , m 

This equation is also RG-invariant. 

By inspection one sees that this equation is satisfied when D =Aj 4(|k|) = -D =j 4a(IpI)> 
which occurs at 

K T = -1, (49) 

or equivalently, at 

k a = (s- 2)/4. (50) 



For this value of k A , the bounds (|44[1 read 

(s - 2)/2 < k v < s/2. (51) 

This expression for k a vanishes for s — 2, and is negative for s < 2, but by the 
inequality (|42|l we have k a > (This expression also violates the horizon condition for 
s < 2.). So the obvious solution to (148 \ , namely kx = —1, does not give a consistent 
solution for s < 2. 
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4.4 Determination of critical exponents 

We substitute the value just obtained, k& = (s — 2)/4, into (|47|) and obtain an equation 
for Ky , 



(i^)r(s/2- Kv )r( Kv - s /2 + i) r(^) 



r(i + Kv ) r( s -K V ) r(^+2)' 

To simplify this equation, we use the identity 

r(s/2 - k v ) r( K y - s /2 + 1) = 



(52) 



sin[7r(s/2 — Ky)] COs(7TKy) 
Here the shifted variable Ky is defined by 

k' v =k v -(s- l)/2, (53) 

and the bounds (|51|) read 

- 1/2 < Ky < 1/2. (54) 
The shifted variable Ky is determined by 

cos(W y ) r[(4i) + «'y] r[(4±) - Ky ] (2^2) r(^) 



6-s^ 



(55) 



The l.h.s. is symmetric about Ky = 0, so if k' v is a solution so is —k' v . Note that Ky = 
is equivalent to Ky = or g 2 DjA A = ( k 2)('s t +i)/2 ■ This describes a linearly rising 
color- Coulomb potential in spatial dimension s. 

Our results are not unreasonable at the physical value s = 3. Indeed, at s = 3, 
the r.h.s. of the last equation has the value ^^^^^ — 21/64 ~ 0.328, while in 
the allowed interval (|54|l . the l.h.s. of (|55|l has a maximum at Ky = where it has 
the value 1/7T ~ 0.318. Thus for s = 3 the two sides agree to within 3%, at Ky — 0, 
which describes a linearly rising color-Coulomb potential, in accordance with numerical 
simulations [30] . 

Note, however, that the maximum of the l.h.s. of (|55p in the allowed interval is 
1/7T, and 1/7T < 21/64. Hence, strictly speaking, there is no real solution at s = 3. As 
s decreases from s = 3, a real solution develops at Ky = and s = s c , where s c is a 
critical dimension that satisfies (|55[1 at Ky — 0, 



r 2 (^±i) _ (^) r(^+2) 



(56) 



namely 

s c « 2.9665. (57) 

Recall that Ky = is equivalent to Ky = (s c — l)/2, and this corresponds to a linearly 
rising color-Coulomb potential at s = s c . As s decreases below the critical dimension 
s c , the solution Ky = bifurcates into two solutions that approach Ky = ±1/2, as s 
approaches s — 2. Although neither of the two branches corresponds to a linearly rising 
color-Coulomb potential, it is intriguing that the average of the two solutions ±k' v is 
Ky — 0, or Ky = (s — l)/2, and this describes a linearly rising potential for all s. 

We have seen that there is no consistent solution for s < 2, and in fact the r.h.s. 
of (|55[) is negative for s < 2, whereas the l.h.s. is positive for k' v in the allowed inter- 
val (]54p . The question, whether vertex corrections could provide a consistent solution 
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will be discussed in section [5] Unfortunately, it turns out that this is not the 
least, at the next-to-trivial level. 

Despite the near success in s — 3 spatial dimensions, we consider the failure at 
s = 2 to indicate a serious flaw of the present truncation scheme. 



4.5 Infrared fixed point of running coupling constant 

In [32], a running coupling constant was defined in gauges that interpolate between 
the Landau gauge and the Coulomb gauge. It was found that in Landau gauge and in 
all interpolating gauges this running coupling has an infrared fixed point 

"interpolating (0) = -jj-. (58) 

In Coulomb gauge this infrared fixed point is given [32] by the formula (in present 
notation) 

4b 2 A b T 

"coulomb (0) = — ■ (59) 

Stv 

However the Coulomb gauge is a singular limit of the interpolating gauge, and the 
methods used in |32j were not sufficient to determine the numerical value of the infrared 
fixed point in Coulomb gauge. However our (near) solution of the Coulomb-gauge DSE 
provides a value. We take spatial dimension s — 3, for which in (|50[) has the value 
KA = 1/4, and substitute these values into (|43|l . This gives 

9 i 47V 

(brblr 1 = ^f, (60) 

and we obtain 

7 

a C oulomb(0) = J7- (61) 

It is not clear what is the origin of the difference between interpolating gauges and the 
Coulomb gauge. There could be a discontinuity but, if so, the discontinuity is not very 
large. Note, however, that different approximations are made in the two calculations. 
Indeed here we used the phase space representation for the Coulomb gauge, whereas in 
interpolating gauges the configuration space representation was used, so the truncation 
schemes are not identical even though both use tree-level vertices. Furthermore, lattice 
investigations of the Coulomb gauge limit of the interpolating gauge suggest also a 
smooth limit, except at zero momentum |20j . 



5 Vertices beyond tree-level 

Since the truncation proposed ultimately has failed, it is necessary to reconsider it. The 
simplest way is to include vertex correction, and check, whether they can be arranged 
such as to arrive at a consistent solutioijj. The most primitive vertex modification 
in this truncation is to keep the tree-level Lorentz structure, but to include a scalar 
dressing function for each vertex appearing. In particular, these corrections will also 
be of instantaneous type. This implies that only the DSEs (|38ll . (|40l) . and (|45|) are 
affected, as all the other consequences of the truncation remain valid. 



5 Note that a similar program has also bee pursued in Landau gauge in 1311 . 
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Except for one important caveat. One of the central elements in the truncation so far 
was the cancellation of the ghost loop and the <^4o-l°°P in the DSE for the transverse 
gluon propagator. That is true at one-loop order in perturbative calculations, since the 
propagators are identical, apart from a factor i. This identity persists beyond one-loop 
perturbation theory also in the instantaneous approximation employed here. However, 
this alone does not guarantee the cancellation of the loops. Either the vertices must also 
coincide, or their differences must arrange, in a subtle way, to allow the cancellation, at 
least, for the purpose at hand, the cancellation of the leading infrared part. To achieve 
self-consistency, the cancellation has therefore to be built in into the ansatze for the 
vertices. 

Therefore, the DSEs to be satisfied are a system of equations for the three unknown 
functions 

D =A t a t Dq Dqq 

and in addition also have to act as constraint equations which restrict the form of 
possible vertex modifications for the three vertices 

pA T cc pA T A <fi pA T A TT T 

The calculations in section Uto obtain the DSEs (|38|l . (|40|l . and ()45|) can still be made 
in the same way, it is only necessary to keep the scalar dressings of the vertices at each 
step. 

To solve the equations, the same ansatze (|39p for the propagators will be made. 
For the scalar vertex dressings the ansatze 

will be made. Herein, anticipating a ghost- ant i- ghost symmetry as in Landau gauge 
[11] , the exponents for the ghost- and anti-ghost-leg of the ghost-gluon vertex have 
been set equal. Furthermore, to guarantee exact cancellation of the ghost diagram 
in the transverse gluon equation by the Aq4> diagram, the parameters of the ghost- 



r AT5c ( P , q ,k) 


= a a p 


r ATA ^( P , q ,k) 


= a a p 




= a (j p 



gluon vertex and the A T Ao<ji have been set equajj. Finally, a A will be set to zero, 
to incorporate the corresponding non-renormalization theorem, which at least exists 
for the limiting interpolating gauge to Coulomb gauge [32]. Of course, the latter is an 
assumption at this stage. The desired solution should have the property ny = 1 in four 
dimensions (s = 3) and Ky = 1/2 in three dimensions (s = 2) to obtain an appropriate 
potential. 

With power-law ansatze and the dimensional renormalization, the integrals can be 
evaluated analytically. The ghost equation (|40p then yields 

^2+2ka _ — 2-2K r +2a T -2K / i x (ftO) 



b 2 A bTa a 

g 2 N(s - 1) T(f - 1 - k t + a T )r(§ - K A )r(2 - f + k t - a T + k a ) 
X 2 s + 1 tt s r(2 + K T -a T )r(s-l-K T + a T -K A )r(l + K A ) 



° If this constraint is not included, it is possible to construct a solution with Ky = 1 with a 
bare ghost-gluon vertex. See also below. 
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where again the horizon condition has been implemented. The equation (|45[) . leads to 

-t. p v ^ = p + p 1 1 x (63) 

b z A b T a a 

g 2 N(s - 1) f(-l + § + q t - K T )r(f - Ky)f(2 - f - a T + n T + «y) 
2S+17J-S f(2 - QT + K T )r(s - 1 + a T - k t - k v )F{1 + k v ) 

Finally, the equation f)38f) yields the lengthy expression 

Ab 2 T p~ 2 ~ iKT - 1 = (64) 

g 2 Na,f 3 b T bo T(f - 1 + /? - b v ) _ 4+s+2/ 3 0+ 2/3 x +2/3 t -kt-kv- 
16tt- 2s r(l - Pq + k v ) P 

((2 2 - s 7r s r(| -l + pr- K T )r{2 - 1 - p - /3r + + k V ) 

x (6 - 3s + 4/3^ + A)(4/3t - 2 - 4k t ) + 4/3 T (-3 + s - 2«* - Ky) + 2«y 
+4k t (3 - s + k t + Ky))) 
1 



+ 



r(2 - pr + K T )r( s - l + p + Pt - «r - k v ) 
r(f + Ar + K T )r(-i -i-Po-p t -k t + k v ) 



4 s r(-/3 7r - Kr)r(i + s + p + P-k + kt - K v) 

(2 2+s n d pi(3, + 2/3* + 2k t )(s + 2/3* + 2k t ) + {2n) s (s + 2/3^ + 2k t ) 

x(-2 + 5S + 4/3 2 - +4k t (1 + s + k t - K V )+4p 7T (l + s + 2k t - k v ) - 6k v ) 
X(2 + s + 2Pk + 2k t - 2k v ) + 4/3 (2 1+s s(l + 2s)n s + 3s(5 + s)(27r) s (y3 w + n T ) 
+2 1+s 7r s ((/3 7r + k t )(2 + 4/3 2 + k t (7 + 5d + 4k t ) + p n (7 + 5s + 8« T )) 
-(3 + 2/3,, + 2k t )(s + 2p n + 2K T )rcy))))) ■ (65) 
The characteristics of these equations are better visible in the form 

7TT P = f(s,a T - K T ,K A )p (66) 

1 2-2k v -+4k z i 2 . /./ \ s— 2+2ar— 2«T_2«y 

P =P + /(s,aT -«r,«v)P (67) 



4&|p- 2 - 4Kr = i + h(s,p T ,p ,Pi r ,K V ,K T ) p - i+s+2f3o+2f3 <' +2 f 3T - KT - K Xe8) 

(69) 

The appearance of the same function / in both equations (|66|l and ()67p is not a co- 
incidence but due to the same type of coupling and integral kernel appearing. This is 
entirely due to the truncation, which implied tree-level tensor structures and the same 
dressing function for both. Of course, an immediate solution would be then Ky = k^, 
but unfortunately, this is forbidden by integral convergence (|42[) . 

Furthermore, the appearance of the dressing of the ghost and ^-vertices only shift 
the gluon exponent in the respective equations by ax, but yields no modification of 
neither the exponent nor the pre- factor consistency conditions. Only in the Ax equation 
(|69|) a consequence would be possible, but it will turn out that this is of minor relevance. 

Let us start with the case s — 3. First of all, the consistency of the exponents on 
both sides of (|67p yields the same relation as from the equation (I66p . 

s — 4 

= — K,x + a T — 2k&, 
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d 

-0.02 r 




0.0 0.1 0.2 0.3 0.4 0.5 




Fig. 4 Top panel: The function d = f(3, -1/2 - 2ka,ka) - /(3, -1/2 - 2kaA) within the 
permitted range < ka < 1/2 for Ky = 1. A zero crossing would indicate a solution to the 
consistency condition. Bottom panel: The function d = /(3, —1/2 — 2re_/\, Ka,) — /(3, —1/2 — 
2ka, Ky). The permitted ranges for and Ky is the lower left triangle of the plot, i. e., the 
part including the origin up to the wall-like structure. 
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as Ky drops out of both sides of l|67[) . For the tree-level term in equation (|67[) to be 
subleading, it is required that 2k ^ < Ky, if both are positive. Trying to solve for the 
pre-factor consistency condition with this limitation, and the limitations < ka < 
min(l,Ky/2) from integral convergence of ([66} and (|67[) and Ky < s/2 from equation 
(|67l) . yield that again no solution in three spatial dimensions exists. This is directly 
visible from the figure [4] left panel. Only if one allows the ghost and the AQ<f> vertices 
to vary independently, is it possible to find a solution. This is then also consistent with 
the equation (I69|) . However, it is no longer self-consistent. So this is the unfortunate 
case, where a solution by power-counting exists, but there is however no solution of 
the system. 

From figure [4] right panel, it is clear that relaxing the restriction Ky = 1 is not 
sufficient to find a solution: In all cases permitted by integral convergence, no solution 
exists. Only in the sense of section l4~4l if one would permit a deviation of s from 3, a 
solution could be found in the same manner. 

In principle, it is possible to relax the integral convergence condition for the equa- 
tion (|67[) . and admit a logarithmic divergence. In this case, Ky = 2k a is also possi- 
ble. This automatically yields that by power-counting all terms scale as the tree-level 
term. It can therefore no longer be dropped. However, a possible renormalization con- 
stant modifies the value of the tree-level constant also by finite parts. This in turn 
modifies the consistency condition even further. It is then possible to find a solution 
for the equations (|66|l and (|67|) . Selecting, e. g., that the finite part of the tree-level 
term is not modified, and thus a MS-like renormalization scheme, yields a solution as 
(k a ,k v ) « (0.66645, 1.3329). 

However, even then the case Ky — 1 remains pathological, because the function / 
diverges for the required value k a = 1/2. 

In addition, the case Ky = 2k & would be due to IR-UV mixing, as the value of 
the exponent becomes entirely determined by the renormalization prescription, which 
seems undesirable. Still, it would be a solution, and the remaining equation, (|69J> , can 
then be solved with the freedom of four additional parameters, Pt, Po, P-k, and kj<, 
since only the combination ax — kt is restricted by (|66|l and (|67|) . 

Hence, there seems to be no pure infrared solution to the system in this truncation 
in a self-consistent manner at s = 3. Introduction of additional vertex tensor structures 
or relaxing some assumptions may be sufficient, but this is currently unclear. 

Finally, note that if one accepted a non-self-consistency in the vertex ansatze, there 
would be solutions with Ky = 1 for appropriate different choices of the exponents in 
(16611 and (|67l) and of the other exponents in (|69l) . In this case, even a bare ghost-gluon 
vertex would be permitted. 

The situation in three dimensions, s — 2, is again different. There is still no solution 
for Ky < 1/2. However, for Ky > 1/2, there is a solution, with k a starting to grow 
from zero at Ky — 1/2. By adjusting ay appropriately, it is then possible to have a 
vanishing gluon propagator even at Ky only marginally greater than 1/2. In addition, 
the equation (|69p is then easily solved, if h is positive for some set of the parameters, 
which can be be achieved even at Pt = /3 T = Pq — 0. However, a massive gluon is then 
impossible, as for kj< — — 1 the function h vanishes. 

So, neither in four nor in three space-time dimensions in this truncation scheme 
does the desired solution exist. In three space dimensions no solution exists at all, and 
in two space dimensions only a solution with Dqq that is too divergent. 
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6 Summary 

In the present work we have studied the infrared limit of the DSEs in the first-order for- 
malism of Coulomb-gauge Yang-Mills theory. Our truncation was based on the assump- 
tion of dominance by the instantaneous contributions. However, it turned out that, even 
taking into account a certain class of non-trivial vertices, no solution was found which 
directly exhibited the desired linear rising Coulomb string tension. However, once this 
condition is relaxed, solutions can be found. At non-vanishing temperatures it is quite 
reasonable to alleviate some of the conditions, corresponding approaches are discussed 
in [33] , 

We also think that the BRST-type on-shell operator r introduced in section [2] may 
be useful in future studies of the Coulomb gauge. 

The results presented here will be useful in constructing other truncation schemes, 
which may be better able to capture the infrared dynamics in this formalism. From 
our discussions it has become clear that such a truncation scheme has to go beyond 
the instantaneous approximation. This then also constitutes the main result of our 
investigation. 
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